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Abstract 



(^ i We introduce a new tool, the Deformed TBA (Deformed Thermo- 

^sj ' dynamic Bethe Ansatz), to analyze the monodromy problem of the cu- 

bic oscillator. The Deformed TBA is a system of five coupled nonlinear 
integral equations, which in a particular case reduces to the Zamolod- 
chikov TBA equation for the 3-state Potts model. Our method gener- 
alizes the Dorey-Tateo analysis of the (monomial) cubic oscillator. We 
(~| . introduce a Y-system corresponding to the Deformed TBA and give it 

m^ ' an elegant geometric interpretation. 

1 Introduction 

en ' 

^ ' This is the first in a series of papers that we dedicate to studying 

>0 , the exact theory of the direct and inverse monodromy problem for the 

'^ ' cubic anharmonic oscillator 



^"^^^^ -V{X;a,b)il;{X) , V{X;a,b) ^ AX^ - aX - b . (1) 



o, 

f^ ^ The main purpose of the present paper is to introduce a novel instru- 

ment of analysis, that we call Deformed Thermodynamic Bethe Ansatz 
(Deformed TBA). 

The monodromy problem of the cubic anharmonic oscillator is a 
. _ fundamental and rather interesting problem in itself. Moreover, it is 

jrt ' deeply interconnected with the study of singularities of solutions of the 

first Painleve equations (see author's papers |Mas| . |MaslO| l. 

Monodromy problems for anharmonic oscillators, and especially 
eigenvalue problems, has been intensively studied since early '70 years 
with a wealth of different methods: see for example the seminal papers 
|BW68| . |Sim70| . |Vor83| . [BBM] . [DT99] . |BLZ01| . [EG09] . 

Computations are often performed by means of the complex WKB 
methods (see for example [DTOOj . |BBM+Ol| and author's papers |Mas) . 
|MaslO| ) or refined asymptotic expansions (see |JSZJ09) and references 
therein). The first breaktrough towards an exact evaluation of the 
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monodromy problem is the work of Dorey and Tateo |DT99| : they an- 
alyze anharmonic oscillators with a monomial potential X^ — E (n not 
necessarily 3) via the Thermodynamic Bethe Ansatz and other nonlin- 
ear integral equations (called sometimes Destri - de Vega equations). 
Subsequently Bazhanov, Lukyanov and Zamolodchikov generalized the 
Dorey- Tateo analysis to monomial potentials with a centrifugal term 
IBLZOlj . 

Here we generalize Dorey and Tateo approach to the general cubic 
potential. In Theorem [5] below, we show that the monodromy prob- 
lem for the general cubic potential is encoded in a nonlinear nonlocal 
Riemann- Hubert problem, which is equivalent (at least for small value 
of the parameter a in ([T])) to the following system of nonlinear integral 
equations that we call Deformed Thermodynamic Bethe Ansatz: 



/-l-oo 
ipiia-a')Aiia'),a,a' eR, ? G Z5 = {-2, . . . , 2} . (2) 

Here 

M^) = E ^'"^Lkia) , Lkia) - In (l + g-^'^t-)) , 



, , 1 V- -^^i^ , ^ yfr 1/3 , ^/3^^ 2/3 ^_,M 
^^, 2*r(ll/6) 4ir(l/6) 

, , VS 2cosh(2cr) , , VS e-i" 



n l + 2cosh(2cr) ' ^^^ ' j^ l + 2cosh(2cr) 

f2{cr) = — — 7—— , (/3-i(cr) = (Pi{-<t) , ^-2{cr) = ^2(-cr) ■ 

TT l + 2cosh(2(T) 

Here the pseudo-energy Sk is related to the logarithm of the k-th 
Stokes multiplier (see Section 3 below). 

For a = 0, see equation ([^^ below, equations ^ reduce to the 
Thermodynamic Bethe Ansatz, introduced by Zamolodchikov [Zam90| 
to describe the thermodynamics of the 3-state Potts model and Lee- 
Yang model. 

The paper is organized as follows. In Section 2 we give a geometric 
construction of the space of monodromy data. We realize the Stokes 
multipliers of the cubic oscillators as natural coordinates on the quo- 
tient W5/PSL{2,C), where W5 is a dense open subset of (P^) (the 
Cartesian product of five copies of P^). Section 3 is devoted to the 
construction of the deformed Y-system. In Section 4 we derive the De- 
formed Thermodynamic Bethe Ansatz. For convenience of the reader, 
we explain the basic theory of cubic oscillators (Stokes sectors, Stokes 
multipliers, subdominant solutions, etc ...) in the Appendix. 
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2 The space of monodromy data 

Usually monodromy data of equation ([1]) are expressed in terms of 
Stokes multipliers, which are defined by means of a special set of so- 
lutions of the equation. In this section, following Nevanlinna |Nev70| 
and author's paper |Mas| . we study the monodromy data from a geo- 
metric (hence invariant) viewpoint. Here and for the rest of the paper 
Z5 = {-2,...,2}. 

Definition 1. Let {(p,x} ^^ '^ basis of solution of (Qp. 
We call 

Wkiip,x)= lim ^TV^^Uoo, fceZg . (3) 

|argA-2^|<f-e 

the k-th asymptotic value. 

We collect the main properties of the asymptotic values in the fol- 
lowing 

Lemma 1. (i) Let ip' = atp + bx and x' = cip + dx', [ ^ ) "= 



Gl{2,C). Th. 



en 



, , awki^,x) + b ... 

Wk[^ ,X) = 7 T-— 7 • (4 

cwk[^,X) +d 

(a) Wk-iiip,x) = Wk+iiip,x) iff o-k{a,b) = . Here (Tk{a,b) is the 
k-th Stokes multipliers. 

(Hi) Wk+i{ip,x) 7^Wk{(p,x) 

Proof See [Mas] . D 

Definition 2. We define 

W5 = {{Z^2, Z^l, Zo, Zl,Z2), Zfc e C U 00, Zfc ^ Zk±l} . 

The group of automorphism of the Riemann sphere, called Mobius 
group or P5L(2,C), has the following natural free action on W5: let 



a 



b 



c 



d 



T == ] e PSL{2, C) then 



„, , ,az^2 + b az2 + b 

T{z-2, • ■ ■ , 22) = ( —7, . . . , -—:) 

CZ-2 + d CZ2 + d 



After formula ^ and Lemma [TJiii) every basis of solution of ([T]) 
determines a point in W5. After the transformation law (|4]). the 
Schrodinger equation ([T]) determines an orbit of the PSL{2,C) action. 
Hence, we define the space of monodromy data as follows. 

Definition 3. We call space 0/ monodromy data the space of the orbits 
of the PSL{2,C) action and denote it V5: 

V5 = W5/PSL{2,C). 

Remark. The author proved in jMa^ that the space of monodromy 
data V5 is the moduli space of solutions of the first Painleve equation. 
V5 is a smooth manifold and A^o.5 C V5 C A^o.5- 

In the rest of the section we construct natural coordinates of V5 
and interpret them in the spirit of cubic oscillator theory. Define 

Rk:W5 -^ C , fc e Z5 , 

Rk{z^2,---,Z2) = (zi+/£,Z-2+fc;^-l + fe,^2+fe) , (5) 

where (a, b; c, d) — i°_^IL _^x is the cross ratio of four point on the 
sphere. 

We collect the main properties of functions Rk in the following 
Lemma, whose easy proof is left to the reader. 

Lemma 2. (i) 

Rk{Z-2, ■ ■ ■ , Z2) ^ (X),\f{z^2,---,Z2)eW5, 

Rkiz-2, . ■ . , Z2) = iff Zk-1 ^ Zk+1 , (6) 

Rkiz-2, ■■■,Z2) = iff Zk-i = Zk+2 or Zk+i = Zk-2 ■ 

(ii) The functions Rk are invariant under the PS L{2^C) action. Hence 
they are well defined on V5; with a small abuse of notation we 
denote Rk also the functions defined on V5 . 

(Hi) They satisfy the following set of quadratic relation^ 

Rk-2Rk+2 = 1 - i?fe , Vfc e Z5 . (7) 

(iv) The pair Rk,Rk+i is a coordinate system of V^ on the open sub- 
set Rk-2 7^ 0. The pair of coordinate systems {Rk, Rk+i) o,nd 
{Rk+2,Rk-2) form an atlas 0/V5. 

Definition 4. We call any cubic polynomial of the form V{X; a, b) = 
4A'^ — aX — b a cubic potential. The above formula identifies the space 
of cubic potentials with C2 3 (a, b). Through the asymptotic values (Q), 
we define the monodromy map 

r : C2 ^ 1/5 ■ 
T{a, b) is the monodromy data of equation (QJi. 



^almost identical to the relations satisfied by the Stokes multipliers (|24p 



Nevanlinna |Nev70| proved that the monodroiny map T is surjec- 
tive. 

With a smaU abuse of notation we define 

Rk{a,b) = Rkor{a,b) . (8) 

The reader can verify that the Stokes muhiphers (defined precisely 
in the Appendix) are (modulo multiplicative constant) the functions 

Rk- 

ak{a,b) = iRk{a,b) . 

We have thus realized the Stokes multipliers geometrically as natural 
coordinates of the monodromy map from the space of cubic oscillators 
to the quotient W5/PSL{2,C)- 

Remark. The same construction presented here holds for anhar- 
monic oscillators with polynomial potentials of any degree. If n is the 
degree, we denote l/„+2 the space of monodromy data and Rj^ the 
natural functions (defined by formula\^. For example V3 is one point 
and V4 is C* . The functions Rj^ satisfy a system of algebraic rela- 
tions similar to (S^, which will be studied in a subsequent paper. The 
reader should notice that only in the cubic case the functions Rj^ 
coincide with the Stokes multipliers. 

3 Y-system 

Here we introduce the Y-system P^ . which is a fundamental step in 
the derivation of the Deformed TBA. 

We begin with an observation, probably due to Sibuya |Sib75) : 

Lemma 3. Let uj = eJ~s~ and Rk be defined as in 0). Then 

Rk.iuj-^a,ujb) = Rk-2ia,b) . (9) 

Proof. Denote (p{X; a, b) a solution of ^ whose Cauchy data do not 
depend on a, b. It is an entire function of three complex variables with 
some remarkable properties. For any fc € Z5 ip{uj''X;uj'^^a,uj^''b) satis- 
fies the same Schrodinger equation ([T]) as ip{X] a, b). Fix ip{X; a, 6), x(A; a, b) 
linearly independent solutions and define 

Wkia,b) = Wk{(p{X;a,b),x{X;a,b)) , 

Wkia,b) = wi{ip{uj'X;uj'^'-a,u}^^b),xi^''X;uj'^'-a,uj^'-b) . 

Then Wk{u!'^''a,Lu^''b) = il!k-i{a,b). Choose / = 2 to obtain the thesis. 

D 

Due to equations ([7]) and relations 0, the holomorphic functions 
Rk{a,b) satisfies the following system of functional equations, first 
studied by Sibuya |Sib75) 

Rkiuj-'^a, ujb)Rk{uja, uj^^b) ^ 1 - Rk{a, b),\/k^'Lr,. (10) 



We have collected all the elements to introduce the important Y- 
functions and Y-system. 
We fix a e C and define 

Yk{^) = -i?o(w" V eS'') , fc e Z5 . (11) 

Sibuya's equation (jlOp is equivalent to the following system of func- 
tional equations, that we call Deformed Y-system: 

Yk-i{d-t^)Yk+i{'d + t^) = l + Yk{d). (12) 

Remark. If a ~ 0, Y^ = YQ,yk and the system 111^) reduces to just 
one equation, called Y-system, which was introduced by Zamolodchikov 
JZam91}/ in relation with the Lee- Yang and 3-state Potts models. Dorey 
and Tateo ]D T99^ studied the Zamolodchikov Y-system in relation with 
the Schrodinger equation with potential V{X; 0, b) = 4A^ — b. 

3.1 Analytic Properties of Y^ 

In the following theorem we summarize the analytic properties of the 
Y-functions. 

Theorem 1. (i) For any a E C and fc £ Z5, Yk is analytic and i^ 

periodic. If a is real then Yk{'&) = K_fe(79), where — stands for 
complex conjugation, 
(ii) For any a € C and fc G Z5, 

Yk{-d) = exp(Ae^ + Sae*-*^) . (13) 



- 1 



= 0{e-^^^), as Red -^ +00 and \Imi}\ < - , 



Here A = VllilZ^ and B = ^^^^''^ . 
23r(ii/6) 4^r(i/6) 

(Hi) For any a G C and any K eM., Yk{d) is bounded on Red < K. 
Ifa^O, lim^^.oo Yk{d) = i^. 

(iv) If e^~^ a is real non negative then Yk{'d) = implies Imd = i^- 

Ifa^O then Yki-&) = -1 implies ?? = ±if 
(v) Fix e > 0. If a is small enough, then for any fc G Z5, Yk{'d) 7^ 
0, -1 for any ■& G \Imd\ <\-e. 

Proof. (i) Trivial. 

(ii) These "WKB-like" estimates can be found in |Mas| Section 4 or 
in [Slb75] . 

(iii) The boundedness follows directly from the fact that Rk{a,b) is 
analytic in 6 = 0. If a = 6 = 0, then for symmetry reasons one 
can choose ip,x such that Wk = e^~s^ . This implies the thesis. 

(iv) The statement is equivalent to Theorem [3] in the Appendix. 

(v) Since Yk depends analytically from the parameter a, it follows 
from (iv). 

D 



4 Deformed TBA 

This section is devoted to the derivation of the Deformed Thermody- 
namic Bethe Ansatz equations Q u- 

In what follows we always make the following 

Assumptions 1. We assume that there exists an e > such that 

(i) every branch of InYfe is holomorphic on |/to'(?| < -I + e, and 
bounded for i? — >■ — oo. And 

(ii) every branch o/ln(l + y-) is holomorphic on J/ttii?! < +e, and 
bounded for "(9 — > — oo. 

From Theorem [IJiii, v) we know that the assumptions are valid if 
a is small enough. 

We define the following bounded analytic functions on the physical 
strip Imd < j 

e^i^) = Inn (79), (14) 

WON ^o^ y¥r(i/3) ^ V3^r(2/3) ._,2^ 

^ ' ^ ' 2tr(ll/6) 4lr(l/6) 

LkW = ln(l + e~^''('')) . 

Here the branches of logarithms are fixed by requiring 

lim 4(cr + ir) = lim Lk{cr + it) = , V |r| < ^ . (15) 

We remark that by Theorem [IJii) , this choice is always possible. We 
denote Sk pseudo-energies in analogy with the undeformed TBA. 

Due to the Y-system (fT^ . the functions 5k satisfy the following 
nonlinear nonlocal Riemann-Hilbert problem 

4-i(i9 - i^) + 4+i(^ + * J) - Sk{^) = Lk{d) , \Imd\ < e . (16) 

Here the boundary conditions are given by asymptotics (|15p . 

The system (|16|) is Z5 invariant. Hence we diagonalize its linear 
part (the left hand side) by taking its discrete Fourier transform (also 
called Wannier transform): 

Xi{^) = Y. ^'"^4(1?) , SkW = i ^ e-^'^xiW , (17) 

A,(^) = ^ e*^Lfe(^) , Lfe(i9) = i^e-^^AKi9). 
fcGZs ;gZ5 

The above defined functions satisfy the following system 

e-'^^Xii^ + * J) + e'^^Xii^ - ^l) - XiW = MW- (18) 

The system of functional equation (jlSp . may be rewritten in the 
convenient form of a system of coupled integral equations ^ . 



^The reader who wants to repeat all the computations below, should remember that 
Er=o'e'^=0. 



Theorem 2. If a is small enough, the functions xi satisfy the De- 
formed Thermodynamic Bethe Ansatz 

(•+00 
Xi(a)=/ ^i{a-a')Ki{a'),a,a' ^n. 



e fined as 


in \lA]n\} and 




\/3 2 cosh(cr) 


¥'o(f ) 


TT l + 2cosh(2cr) 


ViW) 


^/3 e-y 


TT l + 2cosh(2cr) 


</?2(ct) 


V3 e-t- 


TT 1 + 2 C0Sh(2(T) 


<^-i(o-) 


V3 et- 


TT 1 + 2 C0Sh(2(T) 




V3 ei- 


V—2{<^) 


TT l + 2cosh(2(7) 



(19) 



Proof. If a is small enough we know that the Assumptions [T] are valid. 
Hence the thesis follows from system P^ and the technical Lemma 2] 
below. n 

Lemma 4. Let f : \Im§\ < e — > C 6e a bounded analytic function. 
Then for any I G Z5 there exists a unique function F analytic and 
bounded on \Im'd\ < ^ + s, such that 

e~'^F(z9 + i^) + e*^F(i?-i^)-F(i9) = /(i9) , V |/mz9| <e. 

Moreover, F is expressed through the following integral transform 

r+00 
Fi-d+ir) = / ipi{'d+iT-'0')f{d')da' , V|/mi?| < e, |t| < - , (20) 
J-00 ^ 

provided |/?7i(^ + ir — 'd')\ < ^ and the integration path belongs to the 
strip \Imi!)'\ < e. Here ipi is defined by formula Iil9j). 

Proof. Uniqueness: let Fi,F2 be bounded solution of the functional 
equation 

e-^^F,(^+*^)+e'^f,(^-zJ)-^,(^) = f{a) , j = 1,2, \Im^\ < e . 

Their difference G — Fi — F2 satisfies 



e 



~G{i) + i-) + e'^G(i} - i-) - G{§) = . 



Hence G extends to an entire, 2i7r periodic, bounded function. There- 
fore G is a constant and the only constant satisfying the functional 
relation is zero. 



Existence: one notices that ii 6 ^ ±m-|, n G Z then 

e-'^^MS + *|) + e'^^fiie - ^|) - fi{e) = , VI E Z5 . 

Then a rather standard computation shows that the function F defined 
through formula (P(I)l satisfies all the desired properties. D 

Remark. Once the system of integral equations (^ is solved for a G 
M, one can use the same set of integral equations as explicit formulas 
to extend the functions Xi(^) '"^ |/77i0| < ?■. Then one can use the Y- 
system il2\} to extend the Y functions on the entire fundamental strip 
\Im-d\ < ^. 

Remark. While the Y-system equations do not depend on the param- 
eter a (the coefficient of the linear term of the potential AX^ — aX—b), on 
the contrary the Deformed TBA equations depend on it since it enters 
explicitely into the definition of functions A/. 

4.1 The case a = 

If a = then 6k ~ Sq, L^ — L^Vk. Therefore, Sq satisfy the single 
functional equation 

H"^ - ^\) + H^ + ^\) - 'Jo W = ioW , |/mi?| < e . (21) 

Similar reasoning as in Theorem [5] shows that ^o satisfies the follow- 
ing nonlinear integral equation (as it was firstly discovered by Dorey 
and Tateo (DT99]) 

/+00 
ifoia - a') In (1 + cxp (-(<5o(ff') + Ae"' )) , a, a' e R . 
-CXD 

(22) 
r(i/3) 



Here tpo is defined as in formula (|19p and A~^ 



23r(ii/6) 

Equation (P^ is called Thermodynamic Bethe Ansatz and was in- 
troduced by Zamolodchikov |Zam90) to describe the Thermodynamics 
of the 3-state Potts and Lee- Yang models. 

5 Concluding Remarks 

We have given a geometric construction of the space of monodromy data 
of cubic oscillators and we have used such construction to analyze the 
direct monodromy problem. 

We have shown that for small value of the deformation parameter 
(the coefficient of the linear term of the potential) the direct mon- 
odromy problem defines a nonlinear nonlocal Riemann-Hilbert prob- 
lem which is equivalent to a system of nonlinear equations, that we 
call Deformed Thermodynamic Bethe Ansatz (Deformed TBA). Our 
approach can be easily generalized to anharmonic oscillators of any 



order. We will give the details in a subsequent publication (see also 
the Remark at the end of Section 2). 

As it was said in the introduction, this is just the first of a series of 
paper that we dedicate to the Deformed TEA. Our rather ambitious 
purpose is to use the Deformed TEA to construct analytical and nu- 
merical tools to effectively solve the monodromy problem of the general 
cubic oscillator. We plan to pursue our study of the Deformed TEA 
equations in different directions to achieve such a goal. 

In particular, together with A. Moro and R. Tateo we are imple- 
menting an algorithm to solve numerically the equations and to un- 
derstand for which values of the parameter a the Assumptions [1] fail. 
Consequently, we will modify the Riemann-Hilbert problem for taking 
into account the multivaluedness of the functions Ai,xi- 

At the same time, we are going to study analytically the conver- 
gence of successive approximation schemes to solve the Deformed TEA. 

As it is known from author previous works |Mas| . [MaslO| (see also 
|KT05| ). cubic oscillators are deeply interconnected with the first Painleve 
equation and in particular with the distribution of poles of its solu- 
tions. As an application of the theory developed here, we are going to 
investigate the poles of special solution of first Painleve equation. In 
particular we will study the tritronquee solution, since it is relevant 
to describe the singular limit of the focusing Nonlinear Schrodinger 
equation |DGK09| . jBTlOj . 

6 Appendix 

The reader expert in anharmonic oscillators theory will skip this Ap- 
pendix; for her, it will be enough to know that we denote (7k{a, b) the 
k — th Stokes multipliers of equation ^. Here we review briefly the 
standard way, i.e. by means of Stokes multipliers, of introducing the 
monodromy problem for equation ([T]) . All the statements of this sec- 
tion are proved in Appendix A of author's paper |Mas| and in Sibuya's 
book |Sib75| . 

Lemma 5. Fix fc e Z5 = {—2, . . . , 2}, define the branch of A2 by 
requiring 

lim ReA^ = +00 

A-->oo 
arg A=2|i 

while choose the branch of X^ globally on the complex plane minus the 
negative real axis such that it is positive on the positive real axis. Then 
there exists a unique solution yk{X;a,b) of equation {I]j such that 

hm ^^(^i^^^l,Vs>0. (23) 

Definition. We denote {A € C, |arg A - ^| < f } the k-th Stokes sec- 
tor. We denote y^, defined in Lemma above, the k-th subdominant 
solution or the solution subdominant in the k-th sector. 
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From the asymptotics (|23p . it follows that j/^ and yk+i are linearly 
independent and the following equations hold true: 



yk~i{X;a,b) = yk+i{X;a,b) + akia,b)ykiX;a,b) , 



(24) 



-tO-k+3 



1 + (Tfc(a, b)(7k+i{a, 6) , Vfc e Z5 



Definition. The entire functions (Jk{a,b) are called Stokes multipli- 
ers. The quintuplet of Stokes multipliers ak{a,b)^k G Z5 is called the 
monodromy data of equation (Q}). 

6.1 PT-symmetric anharmonic oscillator 

The eigenvalues problem ak = 0,k E Z5 is PT symmetric if w'"'a G M 
(the study of PT symmetric oscillators began in the seminal paper 
HHM]). Dorey, Dunning, Tateo [DDTOlj and Shin |Shi02| proved the 
following result about PT symmetric cubic oscillator. 

Theorem 3. Fix fc G Z5. Suppose (Tk{a,b) = and tJ^^a is real and 
positive. Then ur^^b is real and negative. 
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